Interplay of electron-phonon interaction and strong correlations: DMFT+E study 
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We perform investigation of Hubbard model with interaction between strongly correlated con- 
ducting electrons on a lattice with Debye phonons. To solve the problem generalized dynami- 
cal mean-field DMFT-I-E method is employed with "external" self-energy Ep^ corresponding to 
electron-phonon interaction. We present DMFT-|-Eph results for densities of states and kinks in 
energy dispersions for a variety of model parameters, analyzing the interplay of recently discovered 
kinks of purely electronic nature and usual phonon kinks in the electronic spectrum. 
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I. INTRODUCTION 

The problem of the interplay of strong electronic corre- 
lations with electron-phonon interaction is of central im- 
portance in the physics of highly correlated systems. Ac- 
tually there is rather long history of such studies, e.g. one 
of the most popular models for electron-phonon interac- 
tion (EPI) in strongly correlated systems is the so-called 
Hubbard-Holstein model (HHM). The Hubbard modelJ. 
itself describes local Coulomb interaction of electrons on 
a lattice including e.g. Mott-Hubbard metal-insulator 
transition. On the other hand Holstein model contains 
local linear displacement-to-density interaction of con- 
ducting electrons with local (Einstein) phonon modest. 

Active investigations of the properties of the HHM 
were undertaken in the framework of dynamical mean- 
field theory (DMFT)'^, which is non-perturbative ap- 
proach with respect to interaction parameters of the 
Hubbard model. Among many others one should men- 
tion DMFT solution of HHM for the case where im- 
purity solver used was the numerical renormalization 
group (NRG) (see for review of DMFT(NRG) appli- 
cations Ref. The mapping of HHM to Anderson- 
Holstein impurity was first performed by Hewson and 
Mayer -. It was shown that using NRG one can compute 
in a numerically exact manner total electron-phonon con- 
tribution to the self-energy of the problem, thus making 
solution of the HHM non-perturbative also with respect 
to electron-phonon coupling strength. One should note 
that self-consistent set of DMFT equations is preserved 
in this approach. 

However, up to now there are apparently no stud- 
ies of strongly correlated electrons interacting with De- 
bye phonons. It is even more surprising in view of the 
widely discussed physics of kinks in electronic dispersion 
observed in ARPES experiments 40-70 meV below the 
Fermi level of high-temperature superconductors^, which 
are often attributed to EPF. To our knowledge prob- 
lem of kink formation on electronic dispersion caused by 
EPI in strongly correlated systems was briefly discussed 
within HHM in papers by HagueS. and KoUer et al^. 

In this paper we report DMFT-I-E results for the 
Hubbard model supplemented with Debye phonons, as- 



sumung the validity of Migdal theorem (adiabatic ap- 
proximation). We consider the influence of Debye 
phonons on the weakly and strongly correlated electrons, 
studying electron dispersion and density of states (DOS), 
in particular close to Mott-Hubbard metal insulator tran- 
sition. We analyze in details how EPI affects electronic 
dispersions in correlated metal and discuss the interplay 
of recently discovered kinks of purely electronic nature 
in electronic dispersion^" and usual phonon kinks in the 
electronic spectra. 

The paper is organized as follows. First we introduce in 
Scc.|niDMFT-KS approach to the model at hand. Then 
in Sec, mil calculated results are presented and discussed. 
Suumary and conclusions are given in Sec. IIVI 



II. DMFT+E COMPUTATIONAL DETAILS 

The major assumption of our DMET-I-E approach is 
that the lattice and time Fourier transform of the single- 
particle Green function can be written as: 
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e + M-e(p)-S(£)-Sp(£) 
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where e(p) is the bare electron dispersion, S(£) is the 
local self-energy of DMFT, while Sp(e) is some "ex- 
ternal" (in general case momentum dependent) self- 
energy. Advantage of our generalized approach is the 
additive form of the self-energy (neglect of interference) 
in Eq. |Jj) iiii2,i3 ^ allows one to keep the set of 
self-consistent equations of standart DMFT'^. However 
there are two distinctions. First, on each DMFT itera- 
tions we recalculate corresponding "external" self-energy 
Sp(/i, e, [S(e)]) within some (approximate) scheme, tak- 
ing into account interactions e.g. with collective modes 
(phonons, magnons etc.) or some order parameter fluc- 
tuations. Second, the local Green's function of effective 
impurity problem is defined as 
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at each step of the standard DMFT procedure. 
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FIG. 1: Migdal-like contribution to electron-phonon self- 
energy included into DMFT+Ep/i scheme. 



Eventually, we get the desired Green function in the 
form of (H]), where T,{e) and Sp(e) are those appearing 
at the end of our iteration procedure. 

To treat electron-phonon interaction for strongly cor- 
related system we just introduce Sp(£) = Ttph{s,p) due 
to electron-phonon interaction within the usual Frohlich 
model. To solve single impurity Anderson problem we 
use NRG-. All calculations are done at nearly zero tem- 
perature and at half filling. For "bare" electrons we as- 
sume semielliptic DOS with half-bandwidth D. 

According to the Migdal theorem in adiabatic 
approximatioi*i^ we can restrict ourselves with the sim- 
plest first order contribution to Sp/i(e, p), shown by dia- 
gramm in Fig. [TJ The main advantage of this is possibil- 
ity to neglect any order vertex corrections due electron- 
phonon coupling which are small over adiabatic param- 



eter < 1. 
written as 
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Contribution shown in Fig. [T] can be 



uj"^ - w5(k) + i5 



At - £(p + k) - T.{e + uj)- Typhis + p k) 



(3) 



where g is the usual electron-phonon interaction con- 
stant, wo(k) is phonon dispersion, which in our case is 
taken as in the standard Debye model 



o.o(k) = u|k|, |k| < 



(4) 



Here u is the sound velocity and is Debye frequency. 

Actually Ep;i(e, p) defined by Eq. ([3]) has weak mo- 
mentum dependence which we can omit and continue 
only with significant frequency dependence. For the De- 
bye spectra @ Eq. ^ can be rewritten as (cf. similar 
analysis in Refii^) 
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with 
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where — e— I](e) — S]p^(£) and ujc = ppu is a character- 
istic frequency of the order of ujd ■ For the case of semiel- 
liptic non-interacting DOS No{e) with half-bandwidth D 
we get: 
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FIG. 2: Comparison of DOSes obtained within standard 
DMFT (dashed lines) and BMFT+'Eph (solid lines) methods 
for strong (upper panel, U/2D=1.25) and weak (lower panel, 
U/2D =0.625) Hubbard interaction regimes. Dimensionless 
electron-phonon coupling constant A=0.8. 



It is convenient to introduce the dimensionless electron- 
phonon coupling constant as^^: 



(8) 



To simplify our analysis we shall not perform fully self- 
consistent calculations neglecting phonon renormaliza- 
tion due to EPI^^, assuming that the phonon spectrum 
^ is fixed by the experiment. 



III. RESULTS AND DISCUSSION 

Let us start from comparison between pure DMFT 
and DMFT-fSp/i DOSes for strong (U/2D=1.25) and 
weak (U/2D=0.625) Hubbard interaction presented in 
Fig. [5] on upper and low panels correspondingly. Di- 
mensionless EPI constant ([5]) used in these calculations 
was A=0.8 and Debye frequency W£)=0.125D. In both 
cases we observe some spectral weight redistribution due 
to EPI. For U/2D=1.25 (upper panel of Fig. HD we 
see the well developed three peak structure typical for 
strongly correlated metals. In the energy interval ±ujo 
around the Fermi energy (which is taken as zero en- 
ergy at all figures below) there is almost no difference 
in the DOS quasiparticle peak line shape obtained from 
pure DMFT and DMFT-hSph. However outside this in- 
terval DMFT+'Eph quasiparticle peak becomes signifi- 
cantly broader with spectral weight coming from Hub- 
bard bands. This broadening of DMFT-fSp/i quasipar- 
ticle peak leads as we show below to inhibiting of metal 
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FIG. 3: Sequence of DOSes obtained within standard DMFT 
(dashed lines) and DMFT+Ep/i (solid lines) methods close to 
metal-insulator transition (from top-left to bottom right) with 
A=0.8. 



to insulator transition. In the case of U/2D=0.625 ttiere 
are no clear Hubbard bands formed but only some "side 
wings" are observed. Spectral weight redistribution on 
the lower panel of Fig. [2] is not dramatic, though quali- 
tatively different from the case of U/2D=I.25. Namely, 
main deviations between pure DMFT and DMFT+Sp/j 
happen in the interval iwu, where one can observe kind 
of "cap" in DMFT+Sp/i DOS. Corresponding spectral 
weight goes to the energies around ±U, where Hubbard 
bands are supposed to form. 

In Fig. [3] we compare the behavior of pure DMFT 
and DMFT-I-Ep,, DOSes for different U/2D values 
close to Mott-Hubbard metal-insulator transition. For 
U/2D=1.56 both standard DMFT and DMFT+T,ph pro- 
duce insulating solution. However there is some differ- 
ence between these solutions. The DMFT+Sp/i Hubbard 
bands are lower and broader than DMFT ones because of 
additional interaction (EPI) included. With decrease of 
U for U/2D=1.51 and 1.47 we observe that DMFT-hSp/, 
results correspond to metallic state (with narrow quasi- 
particle peak at the Fermi level), while conventional 
DMFT still produces insulating solution. Only around 
U/2D=1.43 both DMFT and DMFT-^Ep,, results turn 
out to be metallic. Overall DOSes lineshape is the same 
as discussed above. Thus with increase of U finite EPI 
slightly inhibits Mott-Hubbard transition from metallic 
to insulating phase. This result is similar to what was 
observed for the HHM in weak EPI regimei^ii^iil. 

For more deep insight into these results let us ana- 
lyze the structure of corresponding self-energies S(e) and 
Sp/i(e). In Fig. 3] we show both real and imaginary part 
of these self-energies. EPI changes E(e) rather signifi- 



FIG. 4: (Color online) Upper panel — comparison of stan- 
dard DMFT self-energies S(e) (dashed lines) with self- 
energies renormalized by phonons and obtained within the 
DMFT-|-Eph approximation (solid lines). Lower panel — EPI 
self-energies Eph(e). Black lines - real parts, red lines - imag- 
inary parts. A=0.8, U/2D=1.25. 

cantly (see upper panel of Fig. [4]). At the same time in 
±ujn energy interval we find that slopes of real parts of 
both self-energies (which determines quasiparticle weight 
in the Fermi liquid theory) are almost the same, while 
imaginary parts are very close to zero. Thus quasipar- 
ticle peaks should be essentially identical in this region 
as we showed above (Fig. [2]). At energies higher than 
Debye frequencies ReS](e) goes steeper with respect to 
Re(I] -I- Sph), making DMFT qusiparticle peak in DOS 
narrower above lud thus providing faster metal to insu- 
lator transition at A=:0. For the case of U/2D=0.625 
(not shown here) pure DMFT self-energy and those with 
the account of EPI are nearly identical. Corresponding 
Sp/i is very close to that obtained due to phonons only 
and shown on lower panel of Fig. [5] with dashed lines. It 
produces only the "cap" in the DOS around the Fermi 
level mentioned above. One can say also that such a 
"cap" appears in DOS when energy interval 2ujd is much 
smaller than the quasiparticle peak width. 

Now we address the issue of a sudden change of the 
slope of electronic dispersion, the so-called kinks. It 
is well known that interaction of electrons with some 
bosonic mode produces such a kink. In the case of 
EPI typical kink energy is just the Debye frequency ujd- 
Kinks of purely electronic nature were recently reported 
in Ref. 10. 

The energy of purely electronic kink as derived in 
Ref. [13 for semielliptic bare DOS is given by 

iu* = Zfl{V2-1)D, (9) 
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FIG. 5: Schematic picture of pure electronic kink (panel (a)) 
and phonon kink (panel (b)) in electron energy dispersion 
near the Fermi level, - bare energy dispersion with no 
interactions included; £fi - dispersion around the Fermi level 
with electron interaction included; lj* - electronic kink energy; 
ii>D - phonon kink (Debye) energy; Spe and Spph - shifts of 
dispersion due to pure electronic and phonon kinks. 



where D is 

dReT.) I 



the half of bare bandwidth and Zp^ = 
IS Fermi liquid quasiparticle weight, 
is given by the half-width of 
Schematic 



de \e=eF 

The rough estimate of w 
quasiparticle peak of DOS at its half-height 
pictures of kinks of both kinds close to the Fermi level 
are shown in Fig. [5] Electronic kink (on the right side) 
is rather "round" and usually hard to see. This kink 
is formed by the smaller slope connection of two splited 
branches with initial slope (dashed line) at energy ztw*. 
Far away from the Fermi level both of these branches 
return to the initial dispersion. In contrast the phonon 
kink produces rather sharp deviation from the initial dis- 
persion at but outside ztwD energy interval electron 
dispersion quickly returns to the initial one. 

Our calculations clearly demonstrate that electronic 
kinks are hardly observable on the background of phonon 
kinks (as e.g. on upper panel of Fig. |4]), and special care 
should be taken to separate them by rather fine tuning 
of the parameters of our model. To clarify this situation 
we introduce an additional characteristic of the kink — 
the shift of electron dispersion in momentum space 5p 
at kink energy. From simple geometry we estimate for 
phonon kinks 



Spph = — A 



Vp 



(10) 



where vp the bare Fermi velocity and A was defined in 
Eq. ([5]). For electronic kink the similar estimate is 



LU , ZpL 

v*p Zo 



) - ^A. 



(11) 



where Zq is quasiparticle weight in the case of absence 
of electronic kinks (the same as Zcp defined in Ref. [lol ). 
Velocity Vp is the Fermi velocity of initial dispersion, 
but it can not be just a bare one. As was reported in 
Ref. [l3 electronic kinks can be observed only for rather 
strong Hubbard interaction when three peak structure in 
the DOS is well developed and electronic dispersion is 
strongly renormalized by correlation effects. This renor- 
malization is determined by Ag defined in Eq. (jlOp , which 
can be seen as kind of dimensionless interaction constant. 
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FIG. 6: (Color online) Quasiparticle part of DOSes (see Fig. [2] 
upper panel) (red line) and corresponding real part of addi- 
tive self-energy Re(E -I- Eph) with electron-phonon coupling 
switched off (left panel) and switched on (right panel). A=2.0, 
U/2D=1. 



In the case when both slopes on the Fermi level and out of 
ztw* energy interval are equal there will be no electronic 
kink at all. 

Now we can choose parameters of our model to make 
both kinks simultaneously visible. First of all one should 
take care that u>d ^ i^*- For U/2D=1 with U=3.5 eV we 
get LU* '--^O.ID and a reasonable value of Debye frequency 
is ujd '^O.OID. To make phonon kink pronounced at such 
relatively low Debye frequency (cf. Eq. pH]) ) we have to 
increase EPI constant. So we take A=2.0. Corresponding 
quasiparticle peaks of the DOS together with Re(S-|-Sp/i) 
are shown in Fig. [6l at the left panel EPI is switched 
off, while on the right panel it is switched on. We can 
see that 2w* is approximately width of the quasiparticle 
peak of well developed three peak structure (see upper 
panel of Fig. [2]) and energy position of electronic kinks are 
marked by arrows. On the right side of Fig. [51 where EPI 
is present, phonon kinks at dzLUp) are clearly visible and 
well separated in energy from electronic kink position. 

To demonstrate coexistence of both these types of 
kinks we take a look on energy dispersion of simple cubic 
lattice with nearest neighbors transfers only. Most conve- 
nient is high symmetry direction F— (tt, tt, tt) direction^*^. 
In Fig. [7] dispersion along this direction around Fermi 
level is shown. Black line with diamonds is pure DMFT 
electronic spectrum, while red line with circles represent 
the result of DFMT-|-Sp/i calculations. Electronic and 
phonon kinks are marked with arrows. 

Finally we address to the behavior of phonon kinks in 
electronic spectrum as function of Hubbard interaction 
U. As U/2D ratio grows Fermi velocity in Eq. PH)) goes 
down, so that momentum shift of kink position 6p moves 
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IV. CONCLUSION 




High symmetry direction r-(7C,7l,7l) 



FIG. 7: (Color online) Quasiparticle dispersions obtained 
from standard DMFT (black lines with diamonds) and 
DMFT+Eph (red lines with circles) around the Fermi level 
and along the part of high symmetry direction F — (tt, tt, tt). 



away from pp, while kink energy remains at loo- This 
is confirmed by our direct DMFT+Sp^ calculations pro- 
ducing the overall picture of spectrum evolution shown 
in Fig. E 
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This work is a first attempt to analyze strongly cor- 
related electrons, treated within DMFT approach to the 
Hubbard model, interacting with Debye phonons. EPI is 
treated within the simplest (Migdal theorem) approach 
in adiabatic approximation, allowing the neglect of ver- 
tex corrections. DMFT-|-Sp/i approach allows us to use 
the standard momentum space representation for phonon 
self-energy ([3]), while the general structure of DMFT 
equations remains intact. 

Mild EPI leads to rather insignificant changes of elec- 
tron density of states, both in correlated metal and in 
Mott-insulator state, slightly inhibiting metal to insula- 
tor transition with increase of U. 

However, kinks in the electronic dispersion due to EPI 
dominate for the most typical values of the model param- 
eters, making kinks of purely electronic nature, predicted 
in Refii^, hardly observable. Special care (fine tuning) of 
model parameters is needed to separate these anomalies 
in electronic dispersion in strongly correlated systems. 

We have also studied phonon kinks evolution with the 
strength of electronic correlations demonstrating the sig- 
nificant drop in the slope of electronic dispersion close to 
the Fermi level with the growth of Hubbard interaction 
U. 

We believe that these results may be of importance in 
further studies of the evolution of electronic spectra in 
highly corretaed systems, such as e.g copper - oxides. 
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FIG. 8: (Color online) Quasiparticle dispersions around Fermi 
level with phonon kinks obtained from DMFT-|-Eph calcula- 
tions for different interaction strengths U/2D= 0.5, 0.75, 1.0; 
A = 0.8, a;D=0.1D. 
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